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We analyse the process r~ — > VrTr'n allowing for anoma- 
lous weak charged current magnetic and electric dipole mo- 
ment interactions and determine the effects on the differential 
and total decay rates. Using recent experimental data we de- 
termine the following values for the anomalous magnetic and 
electric dipole moment parameters, respectively: 

K = 0.f6±0.08; and \R\ ^ O.SStoll. 
12.60.Cn, 13.35.Dx, 14.60.Fg, 



I. INTRODUCTION 

In the Standard Model the charged current interac- 
tions of the tau lepton are mediated by the W boson 
with a pure V—A coupling. We consider new derivative 
couplings in the Hamiltonian which are parametrised by 
the parameters k and k, the (CP-conserving) magnetic 
and (CP-violating) electric dipole form factors respec- 
tively These are the charged current analogues of 
the weak neutral current dipole moments, measured us- 
ing Z T~^T~ events [^, and the electromagnetic dipole 
moments measured using Z ^ t^t^j events 
The only limits so far obtained for k and k are derived 
from analyses of the partial widths for t~ i~D£i'r, for 
£ = e,Ai|,|,0. 

In this paper we consider for the first time the effects 
of anomalous charged current dipole moments on tau de- 
cays involving hadrons. We analyse the Vt-tt^tt^ 
process which has largest branching fraction of all the tau 
decay modes. This process is particularly topical due to 
a recently reported difference of 2.2a between the mea- 
sured T~ — > z/T-7r~7r° branching fraction and the (lower) 
value predicted using e+e" — > 7r"'"7r~ data in the neigh- 
bourhood of the p meson resonances and the Conserved 
Vector Current (CVC) hypothesis While this could 
be attributed to a fluctuation, we note that non-zero val- 
ues of K and k would also yield a higher measured value 
for BR(r- v^tt-tt°). 

We present predictions for the differential 
Vt-tt^t:^ decay distributions and the partial width, 
T{t~ i^t-tt^tt"), as functions of k and k. The sensi- 
tivity of the differential distributions is analysed for typ- 



ical samples of i^^tt^tt^ decays in e'*"e^ t^t^ 

events which are reconstructed by the LEP and SLC ex- 
periments. The partial width is compared to the ex- 
perimental measurements of BR(t^ i^^tt^tt") to yield 
quantitative constraints on k and k. 

II. PARAMETRISATION OF ANOMALOUS 
COUPLINGS IN r ^ i^^tt 7r° DECAYS 

The matrix element for the decay r~ Urir^Tr'^ is 
given by 

M = J^i?M: (1) 

where Gp = (1.16639 ± 0.00001) x IQ-^GeV-^ is the 
Fermi constant, Kd = (0.9740 ± 0.0010) is the appro- 
priate CKM matrix element ||l^, and and are 
the leptonic and hadronic currents respectively. The ef- 
fects of anomalous weak charged current dipole moment 
couplings at the tv^W vertex are parametrised by aug- 
menting the usual V — A charged current such that 
is given by 

J'^^^x. (^7^(l-7')-^^(«-*«75))u., (2) 

where a^^^ = i/2[7^, 7"^], is the four- momentum trans- 
fer, and rur = (1777.05t{]:^^)MeV (l| is the tau mass. 
The parameters k and k are in general complex but 
henceforth we assume that k is real, as required by CPT 
invariance. The hadronic current is parametrised as 

i^^ = ^/2F(g2)(5,_52)^ (3) 

where qi denote the four-momenta of the two final-state 
pions and F{q^) is a form factor. 

A convenient choice for the kinematic observables, fol- 
lowing Kiihn and Mirkes is: g^, the invariant mass- 
squared of the hadronic system; cos 9, the cosine of 
the angle between the tau spin-vector and the hadronic 
centre-of-mass direction as seen in the tau rest frame; 
and cos /3, the cosine of the angle between the charged 
pion and the axis pointing in the direction of the lab- 
oratory viewed from the hadronic centre-of-mass frame 
(henceforth referred to as the z-axis). 
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A. Differential Decay Rate, 

dr(T" !/r7r~7r'')/dq^dcos6'dcos/3 



B. Dependence of Apparent Polarisation on 
Anomalous Couplings 



After integration over the unobservable neutrino di- 
rection and the azimuthal angle of the charged pion, and 
neglecting the mass difference between the charged and 
neutral pions, the differential decay rate is given by 



dV 



1 



\V^d\'\F{q')\'SEW 



(47r)3 4m3 

X {[fo + ReiK),h + i\K\^ + k^)f2] 
+ Pr [go + Re{K)gi + Klm{K)g2]} 
dq^ d cos (3 d cos 9 



(4) 



where Pt is the tau polarisation and the factor of Sew — 
f.0194 accounts for electroweak corrections to leading 
logarithm The functions fi and gi (i = 0,1,2) are 
given by 
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and ■0 is the angle between the tau direction of flight in 
the hadronic rest frame and the z-axis. At LEP energies 
the following approximation is valid: 



cost/) 
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1 + 77 cos ( 



with 



(13) 
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F{q^) describes the resonant structure of the two-pion 
invariant mass and the model used to describe it is dis- 
cussed in more detail in the following section. 



The usual determination of tau polarisation from en- 
ergy and angular distributions of decay products of the 
tau depends crucially on the assumed V — A structure 
of the charged current to serve as a polarimeter. Addi- 
tional couplings in T — > p;/ decays will produce measured 
values of polarisation which differ from analyses of other 
T decay modes and the predictions from global fits to 
Electroweak parameters in the context of the Standard 
Model. The observed agreement of polarisation measured 
in — > Vt-tt^t^^ decays with other determinations may 
be used to constrain k and k. 

We first integrate the differential width presented 
above with respect to q^ and cos/3. The q^ dependence of 
F must be explicitly considered prior to this integration. 
F{<f') for the r~ VtH^tt^ channel is dominated by the 
p(770) vector meson with a small admixture of p'(1450) 
and a negligible contribution from the p"(1700), as veri- 
fied by the ALEPH experiment |l^ . We work within the 
context of the Kiihn and Santamaria model |l6|| in which 
the p and p' resonances are each described by a Brcit- 
Wigner propagator with an energy-dependent width 11^] 



mz 



(15) 



where: 



r.(<?^) = r° 



g2 \p{ml)^ 



and 



p{s) = l/2y/s - Ami. 



(16) 



where is a constant. The normalisation is fixed by 
chiral symmetry constraints in the limit of soft meson 
momenta, such that the form factor is given by 



F{q 



2, Bp{q^) + PBp,iq^) 
1 + 13 



(17) 



where /3 = —0.145 ||l6[. The differential width, retaining 
only the 9 dependence, is of the form 



dV^ {A + BPr cos ( 



dcos9 



(18) 



The coefficients A and B depend on k and k and are 
given by 

A = f^,+Re{K)fl + {\K\^ + k^)f^ 

B = gl + Re{n)gl + klm{n)gl (19) 

where, for [3 — —0.145, we obtain the following numerical 
values 
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= f518 6 ± 6 4) • 


10-15 GeV 


(20) 


^ A 7f^ 


: (a) 


fl 

■J i 


= (111. 3± 1.9) • 


10-15 GeV 


(21) 


t> 4.3 




J 2 


= (20.74 ± 0.38) 


10^15 GeV 


(22) 
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(221.6 ±2.04) 


10-15 GeV 


(23) 
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= (-37.11 ± .63) 


• 10-15 GeV 


(24) 
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9\ 


(32.73 ±0.50) 


10-15 GeV 


(25) 






Comparison of the results for the measured "apparent po- 
larisation" from r- Vt-t^^t:^ channel with other chan- 
nels or the results of Electroweak fits would permit con- 
straints to be placed on k and k. 
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FIG. 1. The K sensitivity quantity, (T\/iV, as a function of 
Pr for «: ~ 0, for (a) the Vrii'iv'^ decay mode and (b) 

the T~ —* VrUil" decay mode. 



C. Dependence of the Total Width, r(r 
on Anomalous Couplings 



Integration of Eq. |l8| over cos 9 yields the effect of the 
anomalous couplings on the total rate 

r(r- ^ i/^tt-tt") = r" [l + aii?e(K) + asd^P + k^)] , 

(26) 

which naturally is independent of the polarisation (the 
polarisation term is proportional to cos 9 and therefore 
integrates to zero). The parameters ai and 02 are given 
by 

ai = fl/fl^ 0.202; and (27) 

a2 = fl/f^o= 0-037. (28) 

r"(= /q) represents the partial width in the absence of 
anomalous couplings, i.e. k = k = 0. 



III. SENSITIVITY OF THE DIFFERENTIAL 
DECAY RATE TO /t AND k 

The sensitivity of Eq. || to the dipole moment couplings 
was studied for the case of CP-conserving interaction, i.e. 
K real and k = 0. We consider the quantity 

-1/2 

^"/liV (29) 



N ■ 



as a function of the r polarisation, where a is the ex- 
pected error one standard deviation on k, TV is the num- 
ber of r- — > i/t-tt-tt" decays, and dfl is the elemental 
phase space volume. The choice of the quantity a^/N 
simply reflects the 1 / y/N dependence of the statistical 
error a. 

The distribution / is given by Eq. Q for the — > 
j^T-Tr-Tr" channel and Eq. 11 of Rizzo ^ for t~ UtViI' . 
Figure |l] shows a^fN as a function of Pr for the particular 
case K w for (a) the r- — > Vrir^TT 
(b) the r" — > VrVtt' decay mode 
combined) . 



" decay mode, and 
e or fl, not both 



The T- — !■ u-rT^^T^'^ mode is intrinsically more sensitive 
than the leptonic decay modes and is less dependent on 
the T polarisation. In addition, the branching fraction 
for r- Vt-ti^tt^ is larger than each leptonic channel. 

For example, at the Z peak {Pr « —0.15) with a typical 
sample of reconstructed decays for each LEP experiment 
(- 45 000 r- -> i/^tt-tt" decays and - 30 000 r- -> 
VrVit~ decays) the predicted statistical errors are: 



a I 



0.020, and 
0.065, 



(30) 
(31) 



where detector effects are neglected, apart from their in- 
fluence on the reconstruction efficiency which is reflected 
in the number of decays assumed. The corresponding sta- 
tistical error for the combined e and /i channels is ^--^0.046 
which is more than a factor of two less precise than from 
the T~ UrTr~Tr^ channel alone. 

A practical disadvantage of the semileptonic decay is 
the multi-dimensional character of the distribution func- 
tion. In the analysis of the tau polarisation, this problem 
has been overcome using a single "optimal variable" [p^ . 
Although the optimal variable for the tau polarisation 
is not the optimal variable for k (nor for k), we find it 
still provides distinguishing power. We fit hypothetical 
distributions of the optimal variable for simulated sam- 
ples of 45 000 r- — > i'rTT~Tr^ decays each, generated with 
Pr = —0.15 to represent r's produced at the Z peak. 
Typical errors are Up ^ 0.038 which is degraded com- 
pared to the full multi-dimensional fit but is still statisti- 
cally more sensitive than the combined e and fj, channels. 

The apparent disadvantage of the leptonic channels is, 
however, mitigated by the need to know F{q^) for the 
T- — !■ VrT^^T^^ channel which has a non-negligible sys- 
tematic error, as discussed below. 

In this paper we cannot derive constraints on k and 
k from fits to the differential decay distributions due to 
a lack of the necessary experimental information. We 
can, however, determine constraints from the (intrinsi- 
cally less sensitive) measurements of BR(r- — > Vri^^T^^), 
as described below. 
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IV. CONSTRAINTS ON k AND k FROM 

BR(r~ u^n'n") 

We derive quantitative constraints on k and k by con- 
sidering the likelihood for the theoretical prediction for 
BR(t^ — > i/T-7r~7r°) to agree with the experimentally de- 
termined average value of 



BR{t- 



(25.31 ±0.18)%, 



(32) 



as a function of k and k. T'^ of Eq. ^ is determined from 
e+e^ TT+TT^ data using CVC. A combined analysis 
of all data, allowing for radiative corrections and p — uj 
interference, yields the CVC prediction of jnl 



BR{t- 



= (24.52 ±0.33)%, 



(33) 



where the error includes statistical and systematic un- 
certainties and conservatively allows for a possible sys- 
tematic discrepancy of the DM1 data compared to CMD, 
CMD-2, and OLYA. The experimental value of BR(t- 
Vt-tt^-k^) is higher than the CVC prediction by 2.2 stan- 
dard deviations of the combined error. 



We fix ro to the CVC prediction of BR(r- 



-ttO) 



so that only ai and 02 depend on the description of 
hadronic spectral function. This reduces the sensitivity 
of the results to the details of the hadronic modelling. 

We construct likelihoods as a function of n and k con- 
servatively assuming in each case that the other param- 
eter is zero. The errors are propagated numerically ||l9| ] 
taking into account the error on the experimental mea- 
surement of BR(t~ t^^TT^Tr") (Eq. ^2|) , the uncertainty 
on the CVC prediction (Eq. ^3|), a systematic error of 
0.5% for radiative corrections not included in Sew [0i 
and a systematic error of 0.3% for the effect of p—LO inter- 
ference . Figure || shows the likelihood distributions 
for (a) AC and (b) k. The distribution for k has a single 
peak due to the dominance of the term linear in k in 
Eq. |2^. The distribution for k has two symmetric peaks 
due to the lack of a term linear in k in Eq. 
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therefore 

it is more appropriate to constrain the quantity We 
determine 



K = 0.16 ±0.08, and 
0.. 



3+0.25 
'-0.35' 



(34) 
(35) 



where the errors correspond to the 68% confidence level. 
At the 95% confidence level we constrain n and k to the 
ranges: 



0.00 < K < 0.32 ; and 
0.13 < IkI < 1.33. 



(36) 
(37) 



These results are slightly more than two standard devi- 
ations from the SM expectations of zero which, though 
intriguing, cannot be considered statistically compelling 
evidence of new physics. 
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FIG. 2. The likelihood distributions for (a) k, and (b) R. 

The results may be compared to complementary re- 
sults previously obtained from purely leptonic tau de- 
cays [||, which are k = 0.001 ± 0.008 and k = 0.00 ± 0.16 
or -0.014 < K < 0.016 and \k\ < 0.26 at the 95% C.L. [|. 
These are more restrictive than those we obtain from 
T~ — > i/T-7r~7r° decays. This is partly due to the larger un- 
certainties in the theoretical and experimental values of 
the T~ —> t'T-7r~7r° branching fractions. In principle these 
effects could be counteracted by a higher intrinsic sensi- 
tivity of the T~ — > Vt-tt^tt'^ channel due to larger values of 
ai and 02 relative to the leptonic modes. From our calcu- 
lations, however, we see in retrospect that the numerical 
values for ai and 02 are smaller than their leptonic coun- 
terparts (0.5 and 0.1 respectively). Therefore, if only the 
total rate information is used the z/t-tt 
is less sensitive than the leptonic channels, 
to the higher statistical sensitivity of the r" j^t-tt tt 
channel when the differential decay distribution is anal- 
ysed. 



in contrast 




V. SUMMARY 

We present calculations of the differential and total 
decay rates for the process — > z^^tt^tt*', allowing for 
anomalous charged current magnetic and electric dipole 
moments, k and k respectively. This constitutes the first 
such analysis of a hadronic tau decay mode. 
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The analysis of the differential distributions for the 
T~ i/T-7r~7r° decay mode is found to be statistically 
more sensitive than the corresponding analysis of purely 
leptonic modes, t~ — > VTh'ii~ , irrespective of the tau po- 
larisation. The branching fraction, BR(t^ t^^Tr^Tr'^), 
is also sensitive to k, and k although less so than for the 
leptonic branching fractions. 

By comparison of the measured value of BR(t^ — )■ 
Vt-t^^tt^) with the predictions of CVC, we determine k — 
0.16 ±0.08 and |k| = O.SSlJ^j^. which differ from the SM 
expectations by approximately two standard deviations. 

The values for k and k obtained from BR(r^ 
VrVgi~) are consistent with zero. This could mean that 
the measured result for BR(t^ j/t-tt^tt^) and CVC dif- 
fer only due to a fluctuation, or that there is a theoretical 
or experimental uncertainty which is not correctly taken 
into account. The new e^e^ --> 7r"'"7r~ data in the neigh- 
bourhood of the p meson resonances should reduce the 
experimental uncertainty in the CVC prediction by a fac- 
tor of almost two in 1999 ||ll|. Hopefully these data will 
clarify whether this is a statistical or systematic effect 
or the first indication of some new physics which affects 
hadronic tau decays but not purely leptonic decays. 
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